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Abstract. We study how the Einstein relation between spontaneous fluctuations 
and the response to an external perturbation holds in the absence of currents, for the 
comb model and the elastic single-file, which arc examples of systems with subdiffusive 
transport properties. The relevance of noncquilibrium conditions is investigated: when 
a stationary current (in the form of a drift or an energy flux) is present, the Einstein 
relation breaks down, as it is known to happen in systems with standard diffusion. In 
the case of the comb model, a general relation - appeared in the recent literature - 
between response function and an unperturbed suitable correlation function, allows us 
to explain the observed results. This suggests that a relevant ingredient in breaking 
the Einstein formula, for stationary regimes, is not the anomalous diffusion but the 
presence of currents driving the system out of equilibrium. 
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1. Introduction 

In his seminal paper on the Brownian Motion, Einstein, beyond the celebrated relation 
between the diffusion coefficient D and the Avogadro number, found the first example 
of fluctuation-dissipation relation (FDR). In the absence of external forcing one has, for 
large times t — > oo, 

(x(t)) = , (x 2 (t)) ~ 2Dt , (1) 

where x is the position of the Brownian particle and the average is taken over the 
unperturbed dynamic. Once a small constant external force F is applied one has a 
linear drift 

Sx(t) = (x(t)) F - (x(t)) ~ iiFt (2) 

where (■ ■ -)f indicates the average on the perturbed system, and \x is the mobility of 
the colloidal particle. It is remarkable that (x 2 (t)) is proportional to Sx(t) at any time: 

and the Einstein relation (a special case of the fluctuation-dissipation theorem pQ) holds: 
jj, = (3D, with (3 = 1/ksT the inverse temperature and ks the Boltzmann constant. 

On the other hand it is now well established that beyond the standard diffusion, as 
in (JTJ), one can have systems with anomalous diffusion (see for instance [2J [31 [U [5j [6] ) , 
i.e. 

(x 2 {t)) ~ t 2v with i/ ^ 1/2. (4) 

Formally this corresponds to have D = oo if v > 1/2 (superdiffusion) and D — if 
v < 1/2 (subdiffusion). In this letter we will limit the study to the case v < 1/2. It 
is quite natural to wonder if (and how) the FDR changes in the presence of anomalous 
diffusion, i.e. if instead of ([!]), Eq. (j4j) holds. In some systems it has been showed 
that ([3]) holds even in the subdiffusive case. This has been explicitly proved in systems 
described by a fractional- Fokker-Planck equation [7] , see also [H [9] . In addition there is 
clear analytical [10] and numerical [H] evidences that ([3]) is valid for the elastic single 
file, i.e. a gas of hard rods on a ring with elastic collisions, driven by an external 
thermostat, which exhibits subdiffusive behavior, (x 2 ) ~ t 1 ^ 2 [T2] . 

The aim of this paper is to discuss the validity of the fluctuation-dissipation relation 
in the form ([3]) for systems with anomalous diffusion which are not fully described by 
a fractional Fokker-Planck equation. In particular we will investigate the relevance of 
the anomalous diffusion, the presence of non equilibrium conditions and the (possible) 
role of finite size. Since we are also interested in the study of transient regimes, we 
will consider models with microscopic dynamics described in terms of transition rates 
or microscopic interactions. 

First, we focus on the study of a particle moving on a "finite comb" lattice with teeth 
of size L [13]. In the limit L = oo an anomalous subdiffusive behavior, (x 2 ) ~ t 1 ^ 2 , holds 
and the system can be mapped, for large times, onto a continuous time random walk |13j . 



On anomalous diffusion and the out of equilibrium response function. 



3 



For finite L the subdiffusion is only transient and at very large time t > t*(L) ~ L 2 
one has a standard diffusion: (x 2 ) ~ t. We will see that Eq. ([3]), where in this case the 
perturbed average is obtained with unbalanced transition rates driving the particle along 
the backbone of the comb, holds both for t > t*(L) and t < t*(L) with the same constant. 
This in spite of the fact that the probability densities P(x, t) in the two regimes are very 
different. The scenario changes in the presence of "non equilibrium" conditions, i.e. with 
a drift, which induces a current, in the unperturbed state: the relation ([3]) does not hold 
anymore. On the other hand, in this case it is possible to use a generalized fluctuation- 
dissipation relation, derived by Lippiello et al. in [T3], which gives the response function 
in terms of unperturbed correlation functions and is an example of non equilibrium FDR 
valid under rather general conditions [151 HS1 EH EH EH E3 EQl IZE] . A generalization of 
the Einstein formula was also proved in the framework of continuous time random walks 
in [22] • So we can say that the Einstein relation ((3]) also holds in cases with anomalous 
diffusion when no current is present, but it is necessary to introduce suitable corrections 
when a perturbation is applied to a system with non zero drift. 

In addition we compare the results found in comb models, with those obtained for 
single-file diffusion with a finite number of particles. There we will also consider a non 
equilibrium case, with the introduction of inelastic collisions which induce an energy flux 
crossing the system. Our results suggest that the presence of non equilibrium currents 
plays a relevant role in modifying Eq. ([3]) in stationary states. 



2. Comb: diffusion and response function 

The comb lattice is a discrete structure consisting of an infinite linear chain (backbone), 
the sites of which are connected with other linear chains (teeth) of length L [13]. We 
denote by x £ (— oo, oo) the position of the particle performing the random walk along 
the backbone and with y £ [— L, L] that along a tooth. The transition probabilities 
from (x,y) to (x',y') are: 

W d [(x, 0) ->■ (x ± 1, 0)] = 1/4 ± d 
W d [(x,0) -> (x,±l)] =1/4 

W d [(x,y)^(x,y±l)] = l/2 for y £ 0, ±L. (5) 

On the boundaries of each tooth, y = ±L, the particle is reflected with probability 
1. The case L = oo is obtained in numerical simulations by letting the y coordinate 
increase without boundaries. Here we consider a discrete time process and, of course, 
the normalization Yl(x' y') W d [(x, y) — > (x',y')] = 1 holds. The parameter d £ [0,1/4] 
allows us to consider also the case where a constant external field is applied along the 
x axis, producing a non zero drift of the particle. A state with a non zero drift can be 
considered as a perturbed state (in that case we denote the perturbing field by e), or 
it can be itself the starting state where a further perturbation can be added changing 
d — >■ d+ e. 
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Figure 1. Left panel: (x 2 (t))o/L vs t/L 2 is plotted for several values of L in the comb 
model. Right panel: (x 2 (t)) and the response function 6x(t) for L = 512. In the inset 
the parametric plot Sx(t) vs (x 2 (t))o is shown. 



Let us start by considering the case d = 0. For finite teeth length L < oo, we have 
numerical evidence of a dynamical crossover from a subdiffusive to a simple diffusive 
asymptotic behaviour (see Fig. [TJ) 

/ S Ctl/2 t<t*(L) 

{X W)o " { 2D(L)t t > r(L), (6) 

where C is a constant and D(L) is an effective diffusion coefficient depending on L. The 
symbol (. . .)o denotes an average over different realizations of the dynamics (jSJ) with 
d = and initial condition x(0) = y(0) = 0. We find t*(L) ~ L 2 and D(L) ~ 1/L and 
in the left panel of Fig. [T]we plot {x 2 (t))o/L as function of t/L 2 for several values of L, 
showing an excellent data collapse. 

In the limit of infinite teeth, L — > oo, D — > and t* — > oo and the system shows a 
pure subdiffusive behaviour [23] 

(x 2 (t)) ~ t 1/2 . (7) 

In this case, the probability distribution function behaves as 

( M \ 4/3 

PoMl^'VlW , (8) 

where c is a constant, in agreement with an argument a la Flory [2]. The behaviour (JSJ) 
also holds in the case of finite L, provided that t < t* . For larger times a simple 
Gaussian distribution is observed. Note that, in general, the scaling exponent u, in 
this case v — 1/4, does not determine univocally the shape of the pdf. Indeed, for the 
single-file model, discussed below, we have the same v but the pdf is Gaussian [24J. 

In the comb model with infinite teeth, the FDR in its standard form is fulfilled, 
namely if we apply a constant perturbation e pulling the particles along the 1-d lattice 
one has numerical evidence that 

(x 2 (t)} ~ CSx{t) ~ t 1 ' 2 . (9) 
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In the following section we derive this result from a generalized FDR. Moreover, the 
proportionality between (x 2 (t)) and 5x(t) is fulfilled also with L < oo, where both 
the mean square displacement (m.s.d.) and the drift with an applied force exhibit the 
same crossover from subdiffusive, ~ t 1 / 2 , to diffusive, ~ t (see Fig. (TJ right panel). 
Therefore what we can say is that the FDR is somehow "blind" to the dynamical 
crossover experienced by the system. When the perturbation is applied to a state 
without any current, the proportionality between response and correlation holds despite 
anomalous transport phenomena. 

Our aim here is to show that, differently from what depicted above about the zero 
current situation, within a state with a non zero drift [25] the emergence of a dynamical 
crossover is connected to the breaking of the FDR. Indeed, the m.s.d. in the presence 
of a non zero current, even with L = oo, shows a dynamical crossover 

(x 2 (t)) d ^at 1 / 2 + bf (10) 

where a and b are two constants, whereas 

te d {t) ~ t l '\ (11) 

with Sxd(t) = (x(t))d+ £ — ( x (t))d' at large times the Einstein relation breaks down (see 
Fig. [2]). The proportionality between response and fluctuations cannot be recovered by 
simply replacing (x 2 (t)) d with (x 2 (t))d — {x{t)) 2 d1 as it happens for Gaussian processes 
(see discussion below), namely we find numerically 

([x(t) - (x(t)) d ] 2 ) d ~ a' t 1 ' 2 + b' t, (12) 

where a' and b' are two constants, as reported in Fig. 



3. Comb: application of a generalized FDR 

The discussion of the previous section shows that the first moment of the probability 
distribution function with drift Pd{x,t) and the second moment of Po(x,t) are always 
proportional. Note that in the presence of a drift the pdf is strongly asymmetric with 
respect to the mean value, as shown in Fig. [3] for a system with L = oo. Differently, the 
first moment of Pd+ E (x,t) is not proportional to the second moment of Pd{x,t), namely 
{x(t))d+ £ (x 2 (t)) d — (x(t)) d . In order to find out a relation between such quantities, 
we need a generalized fluctuation-dissipation relation. 

According to the definition (jSJ), one has for the backbone 

W d+S [(x,y) (x',y')} = W d [(x,y) (x',y')} (l + ) ~ W d e^'~ x \ 

(13) 

where W° = 1/4, and the last expression holds under the condition d/W° 1. 
Regarding the above expression as a local detailed balance condition for our Markov 
process we can rewrite it, for (x,y) ^ (x',y'), as 

W**[(x, y) (x', y')} = W d [(x, y) -+ (x', y')\e^'-*\ (14) 
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Figure 2. Response function (black line), m.s.d. (red dotted line) and second 
cumulant (black dotted line) measured in the the comb model with L = oo, field 
d = 0.01 and perturbation e = 0.002. The correlation with activity (green dotted line) 
yields the right correction to recover the full response function (blue dotted line), in 
agreement with the FDR (|15|) . 



where h(e) = 2e/W°. For general models where the perturbation enters the transition 
probabilities according to Eq. ( fl4l) . the following formula for the integrated linear 
response function has been derived [TU [19J [21] 

j|y = m))d X[ £) m))d = \ imWQU - (O(t)x(0)) d - (0(t)A(t, 0)) d ] , (15) 
where O is a generic observable, and A(t,0) = Ylt'=oB(t'), with 

B[(x, y)] =J2( X '~ xW d [(x, y) W, y% (16) 

(x>,y>) 

The above observable yields an effective measure of the propensity of the system to leave 
a certain state (x, y) and, in some contexts, it is referred to as activity [26] . Recalling 
the definitions (EJ), from the above equation we have B[(x, y)] = 2dS yj0 and therefore the 
sum on B has an intuitive meaning: it counts the time spent by the particle on the x 
axis. The results described in the previous section can be then read in the light of the 
fluctuation-dissipation relation ( 1T5|) : 

i) Putting 0{t) = x(t), in the case without drift, i.e. d = 0, one has B = and, 
recalling the choice of the initial condition x(0) = 0, 

Sx {x(t)) e - (x(t)) 1, 2 



m h{e) 2<*« V < 17) 

This explains the observed behaviour even in the anomalous regime and predicts the 
correct proportionality factor, Sx(t) = e/W°(x 2 (t)) . 

ii) Putting 0(t) = x(t), in the case with d ^ 0, one has 

5x d 1 



h(e) 2 



[(x 2 (t)) d -(x(t)A(t,0)) d ]. (18) 
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This explains the observed behaviours fllOp and (11 II) : the leading behavior at large 
times of {x 2 (t)) d ~ t, turns out to be exactly canceled by the term (x(t)A(t, 0)) d , so 
that the relation between response and unperturbed correlation functions is recovered 
(see Fig. |2I). 

hi) As discussed above, it is not enough to substitute (x 2 (t)) d with (x 2 (t)) d — (x{t)) 2 d 
to recover the proportionality with Sx d (t) when the process is not Gaussian. This can be 
explained in the following manner. By making use of the second order out of equilibrium 
FDR derived by Lippiello et al. in [27J EEJ EH], which is needed due to the vanishing of 
the first order term for symmetry, we can explicitly evaluate 



(x 2 (t)) d =(x 2 (t)) + h 2 (d)- 

where A^(t,0) = £^, =0 £ (2) (i') with b(2) = 
— l/25 y>0 . Then, recalling Eq. (IT7|) . we obtain 

1 



:(t)) + hx 2 (t)AW(t,0))t 



(19) 



ZAx' ~ x) 2 W[(x,y) (x>,y>)] 



(x 2 (t)) d -(x(t)) 2 d =(x 2 (t)) + h 2 (d) 



[x\t)) { 



\{x 2 {t)A^ 2 \tM»-\(At)) 



• (20) 



Numerical simulations show that the term in the square brackets grows like t yielding 
a scaling behaviour with time consistent with Eq. ffl~2]) . On the other hand, in the case 
of the simple random walk, one has 5 (2) = -1 and A^ 2 \t, 0) = -t and then 

1 



(x 2 (t)} d -(x(t)}j=(x 2 (t)} + h 2 (d) 



h x \t)) -lt(x 2 (t)) 



4 



[x 2 (t))l 



.(21) 



Since in the Gaussian case (x 4 (t)) = 3(x 2 {t))l and (x 2 (t)) = t, the term in the square 
brackets vanishes identically and that explains why, in the presence of a drift, the second 
cumulant grows exactly as the second moment with no drift. 



0.01 - 



X 




300 



Figure 3. Pd{x, t) in the comb model with L = oo and d = 0.01 at different times. 
Notice that the mean value increases with time mostly due to the spreading, while the 
most probable value remains always close to zero. 
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4. Conclusions and perspectives 

In order to evaluate the generality of the above results, let us conclude by discussing 
another system. Indeed, subdiffusion is present in many different problems where 
geometrical constraints play a central role. In this framework, a well studied 
phenomenon is the so-called single-file diffusion. Namely, we have N Brownian rods 
on a ring of length L interacting with elastic collisions and coupled with a thermal bath. 
The equation of motion for the single particle velocity between collisions is 



where m is the mass, 7 is the friction coefficient, and 77 is a white noise with variance 
(r](t)r](t')) = 2T~f8(t — t'). The combined effect of collisions, noise and geometry (since 
the system is one-dimensional the particles cannot overcome each other) produces a 
non-trivial behaviour. In the thermodynamic limit, i.e. L,N — > 00 with N/L — y p, a 
subdiffusive behaviour occurs [T2j . 

Analogously to the comb model, the case of N and L finite presents some interesting 
aspects. In order to avoid trivial results due to the periodic boundary conditions on 
the ring, it is suitable to define the position of a tagged particle as s(t) = f Q v(t')dt', 
where v(t) is its velocity. For the m.s.d. (s 2 (t)), averaged over the thermalized initial 
conditions and over the noise, we find, after a transient ballistic behaviour for short 
times, a dynamical crossover between two different regimes: 



where a is the length of the rods and D is the diffusion coefficient of the single Brown- 
ian particle [12] ■ Note that the asymptotic behaviour is completely determined by the 
motion of the center of mass, which is not affected by the collisions and simply diffuses. 
Moreover, as evident from numerical simulations, r* ~ N 2 and in the limit of infinite 
number of particles the behaviour becomes subdiffusive, in perfect analogy with what 
observed for the comb model, where the role of L is here played by N. The main differ- 
ence is that, in this case, the probability distribution is Gaussian in both regimes. As a 
consequence of the Gaussian nature of the problem, applying a perturbation as a small 
force F in Eq. ( J22l) . one finds that the Einstein relation is always fulfilled [TT1 ITOl [T8l 130] . 
also for finite N and L (see Fig. H]). Strong violations of the Einstein relation, can be 
obtained, in dense cases, when the collisions between the rods are inelastic so that a 
homogeneous energy current crosses the system [TT] . 

In this note we have considered systems with subdiffusive behaviour, showing that 
the proportionality between response function and correlation breaks down when "non 
equilibrium" conditions are introduced. In the comb model, non equilibrium effects are 
induced by unbalanced transition probabilities driving the particle along the backbone, 
while the single-file model is driven away from equilibrium by inelastic collisions. In 
the first case, the generalized FDR of Eq. (115j) . developed in the framework of aging 



mv(t) 



1 v{t)+r ] {t), 



(22) 




t < T*(N) 
t > T*(N), 



(23) 



On anomalous diffusion and the out of equilibrium response function. 



9 




t 

Figure 4. (s 2 (t)) and the response function Ss(t) for the single-file model with 
parameters: N = 10, L = 10, a = 0.1, m = 1, 7 = 2, T = 1 and perturbation 
F = 0.1. In the inset the parametric plot <5s(t) vs (s 2 (t)) is shown. 

systems [2], can be explicitly written, providing the off equilibrium corrections to the 
Einstein relation. In the second case, the transition rates are not known and another 
formalism must be exploited [18], which requires the knowledge of the probability 
distribution for the relevant dynamical variables of the model. For instance, following 
the ideas of [IT], a distribution which couples the velocities of neighbouring particles 
could be a reasonable guess. Still, the indentification of the relevant variables and their 
coupling in the single-file and other granular systems is a central issue, requiring further 
investigations. 
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